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Abstract 

> ■ 

Measuring the photon energy spectrum in radiative B decays provides essential help 
for gaining theoretical control over semileptonic B transitions. The hadronic recoil 
mass distribution in B^>X u lv promises the best environment for determining \ V u b\. 
The theoretical uncertainties are largest in the domain of low values of the lepton 
pair mass q 2 . Universality relations allow to describe this domain reliably in terms 
of the photon spectrum in B— >X s + , y. A method is proposed to incorporate l/ra& 
corrections into this relation. The \ow-E 1 tail in radiative decays is important in 
the context of extracting \V u b\. We argue that CLEO's recent fit to the spectrum 
underestimates the fraction of the photon spectrum below 2 GeV. Potentially sig- 
nificant uncertainties enter in the theoretical evaluation of the integrated end-point 
lepton spectrum or the B^X u £u width with a too high value of the lower cut on 
q 2 in alternative approaches to \V u b\. 
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1 Introduction 

Since B — > X+y decays represent a loop effect within the Standard Model (SM), they 
are viewed as a window onto New Physics. There is a second motivation that has 
attracted considerable attention: measuring the inclusive photon spectrum allows 
us to infer the motion of the heavy quark inside the B meson. This yields new 
insights into the inner workings of QCD which significantly affect, for example, the 
lepton spectrum in b — > u iv decays. The moments of the heavy quark distribution 
function allow in principle to determine experimentally a number of nonperturbative 
parameters intrinsic to the heavy quark expansion the knowledge of which is of 
primary practical importance |l|, ^ . 

One should note, however, that we are not in the dark about the values of 
these quantities: we know the running beauty quark mass m h and the kinetic ^ 
and chromomagnetic fj, G expectation values with good accuracy due to intensive 
theoretical investigations over the last few years (for a review see Ref. ||). The 
most accurate value for mj has been determined from e + e~ annihilation to beauty 
hadrons close to threshold. Heavy quark sum rules impose tight constraints; in 
particular they ensure /i^. > fi G , which in practice allows for only a limited range 
for the former. Extracting the values of m,b, etc. from the moments or shape of 
energy spectra thus serves as a cross check - and a highly valued one at that - of 
our theoretical control. In lucky cases this can allow to narrow down the theoretical 
uncertainties by excluding some of the corners in parameter space. This comparison 
provides us also with information about the scale of higher order effects which will 
be very important in assessing the numerical accuracy of concrete applications of 
the heavy quark expansion 0]. 

It has always been tempting to extract \V u b\ from the lepton energy endpoint 
spectrum in B — > X u Iv using information on the b quark distribution function in- 
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ferred from the photon spectrum in B — > X s +^f. We will illustrate that theoretical 
problems arise in such a procedure which may severely limit its accuracy. It is more 
profitable to use the information from the photon spectrum to control possible theo- 
retical uncertainties in the hadronic recoil mass spectrum in B — > X u Iv originating 
from low-g 2 kinematics. 

Very recently the CLEO collaboration has presented a new measurement of the 
branching ratio and photon spectrum for B—^Xg+j The bulk of their spectrum 
shown in Fig. 1 is quite similar to the theoretical predictions (Fig. 3 of Ref. ||) made 
already in 1995. Since then our knowledge of m\, and /i 2 has become more precise; 
the value of the latter, which determines the width of the energy distribution is now 
expected in the lower half of the interval used in 1995. This brings the main part of 
the spectrum into even better agreement with the data. One should keep in mind, 
though, that the experimental photon spectrum is additionally Doppler-smeared by 
the nonvanishing velocity /3~0.063 of B mesons in decays of T(45). 



The model spectrum used by CLEO (hereafter referred to as CLEO's fit) assigns 
only a tiny fraction of the photon spectrum to fall below 2 GeV. This would be a 
welcome feature since the background becomes much more severe there; it means 
that the branching ratio for radiative decays can therefore be measured with good 
accuracy. The total radiative branching ratio is little affected by the possible un- 
certainties in this fraction. Yet the photon spectrum below 2 GeV has considerable 
intrinsic interest. Since it is under better theoretical control as far as perturbative 
and nonperturbative effects are concerned, it provides a testbed for studying the 
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Figure 1: CLEO's photon spectrum and the model fit to the data from Ref. |J. 
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onset of local parton-hadron duality. Its possible violation, although of no practi- 
cal concern for integrated widths, can significantly affect the decay rates integrated 
over limited kinematic domains, or interfere with extracting parameters describing 
strong dynamics from higher moments. Besides, data from this kinematical domain 
are important for establishing precision control over the low-g 2 region in B — > X u £u, 
referred to above. 

It is therefore important to have accurate data on the photon spectrum below 
2 GeV despite the experimental challenge it poses. CLEO has performed a prelimi- 
nary measurement of that domain. However, we think that the experimental fit in 
Ref. || significantly underestimates the spectrum that has to exist there. The the- 
oretical framework itself used to translate the observed spectrum into heavy quark 
parameters assumes that a more significant fraction of the photon spectrum falls 
below 2 GeV. We make the case for future data with even higher statistics to clarify 
this issue, and suggest they will show that more than 10% of the photon spectrum 
resides below 2 GeV. An intriguing lesson can, however, be drawn if CLEO's fit 
were to be confirmed. 

The remainder of the paper is organized as follows: after briefly recalling in 
Sect. how the photon spectrum in radiative B decays is obtained theoretically, 
we discuss more specifically the perturbative spectrum mainly responsible for the 
lower tail in Sect. The complete spectrum is considered in Sect. \L% In Sect. |3| 
arguments are given calling for caution when connecting the lepton energy endpoint 
spectrum in B — > X u £u with the photon spectrum in B — > X s + r y, and a simplified 
method to account for a part of subleading corrections is suggested. We then describe 
a way to further reduce the theoretical uncertainties in extracting \V u b\ from the 
hadronic recoil mass spectrum in B — > X u iv in Sect. |], before commenting on some 
issues related to extracting \V u b\ and presenting our conclusions in Sect. ||. 

2 Theoretical expectation for the photon spec- 
trum 

While the integrated width for the inclusive transition B + 7 can be affected 
by the intervention of New Physics, the shape of the prompt photon spectrum is 
more robust against it, at least within non-exotic scenarios. Two classes of effects 
govern the photon spectrum: the lowest order process b — > q + 7 would produce 
a monoenergetic line = b 2m ^ 9 ~ rrib/2 characteristic of two-body kinematics; 
gluon bremsstrahlung adds a continuous hadronic (partonic) mass spectrum all the 
way up to nib thus producing a radiative tail covering the whole range of smaller 
photon energies. 

Gluon field modes with momenta of order Aqcd introduce additional nonper- 
turbative effects: most notably they induce a primordial energy-momentum distri- 
bution of the b quark inside the B meson - an effect that had been introduced on 
phenomenological grounds more than 20 years ago [8]]. The emergence of this 
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'Fermi motion' in the OPE-based QCD heavy quark expansion was pointed out in 
Ref. [0. The related Doppler-smearing further broadens the photon spectrum and 
populates also the gap ^ < E y < ^f- distinguishing quark-level and hadronic kine- 
matics. To leading order in 1/mb the effect of the primordial distribution is given 
by the convolution of the short-distance-generated parton spectrum dr pcrt /di? with 
the heavy quark distribution function F: 

The natural variable for F is the light-cone momentum k + = k Q — k z carried by the 
heavy quark in the parton description. In analogy with the deep inelastic scattering 
structure functions it is often traded for a dimensionless variable x — it with 

A 

A = lim mt ^ 00 (Ms — mb). The function F(x) has support from — oo to 1, and its 
moments are given by the expectation values of local heavy quark operators: 

1 f i 
F(x) dx = 1 , / x F(x) dx = 



oo J —oo 



1 x 2 Fix) dx = , f 1 x 3 Fix) dx = -X.^, (2) 

etc. ifi 2 , p 3 D denote the expectation values of the kinetic energy operator and the 
Darwin term, respectively). At large negative x this function must rapidly vanish. 
In actual QCD the moments and likewise the distribution function itself depend 
on the normalization point, and this dependence is actually far more pronounced 
than in the usual leading-twist functions of light hadrons. The exact separation 
between quantum modes that are short- distance and belong to the hard "pertur- 
bative" kernel, and those which are considered soft and thus incorporated into the 
nonperturbative 'primordial' effects, is to some extent arbitrary and is set by fixing 
a normalization scale fi. In this way both elements on the right hand side of Eq. ([!]) 
are /x-dependent, while the physical spectrum on the left hand side is not. 

Even leaving aside effects of higher order in 1/ m^, predicting the photon spectrum 
would require detailed knowledge of the function Fix), which is not quite realistic. 
Predictions || were obtained assuming a reasonable functional form of Fix) and 
fitting its parameters to reproduce the known and estimated B expectation values 
in Eq. (0). Such a choice for Fix) is not unique, and one might be concerned that 
an essential uncertainty is thus introduced into the spectrum. It turns out, however, 
that its gross features are relatively stable as long as the first three or four moments 
are fixed. 

This stability applies in particular to the lower part of the spectrum, especially 
when one considers the fraction of events with photon energies below a value E: 

*^ = t { bJx, + 1 ) C ae ^- < 3 > 
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The fraction of decay events with E 1 < 2 GeV theoretically comes out about 12%, 
significantly larger than the 5% CLEO's fit would yield and above 8% one would 
obtain literally using the central data points. Yet, it is qualitatively evident that suf- 
ficiently low parts of the spectrum are shaped by emission of relatively hard gluons 
where the perturbative description is adequate and the intervention of nonpertur- 
bative dynamics is insignificant. This was manifest in the analysis of Ref. ||. 

Of course, the uncertainty in the distribution function F(x) can to some extent 
affect the fraction <3> 7 (2GeV) if, say, F(x) has a long tail towards large negative x. 
It is usually missed, however, that such a scenario would imply a rather large value 
for the Darwin operator p|, (and other higher-order ones) by virtue of their relation 
to the moments in Eq. @. This would not go unnoticed in other applications of 
the heavy quark expansion. Yet since our direct knowledge of such expectation 
values is at present still limited, we will not draw conclusions here based on this 
observation. Instead we merely want to point out that a broader distribution F(x) 
yields additional smearing, which tends to increase the low-E^ fraction of decay 
events. To quantify these intuitive assertions we consider in more detail the purely 
perturbative component of the spectrum . 

2.1 The perturbative spectrum 

Throughout this paper we assume the radiatively induced decay to be generated by 
a single local operator 

£wcak = /^<V(1 + 7 5 )&F^ . (4) 

Normalized at the scale mb, this term indeed yields the dominant contribution to the 
b — > s + 7 decay rate in the Standard Model. Other contributions to the transition 
amplitude generated at scales below give sizeable corrections, in particular the 
usual four-fermion operator sb cc with the low-scale annihilation of the cc pair into 
photon and hadrons. However, this mostly modifies the total rate, whereas the 
decay fraction we are interested in is hardly changed. 

Perturbative gluon effects in the b —>■ s-py spectrum contain some theoretical com- 
plexities; for doubly logarithmic Sudakov corrections arise due to the emission prob- 
ability oc ^ where u and k± are the gluon energy and transverse momentum, 
respectively. The gluon coupling grows at small kj_, which on one hand enhances 
corrections and on the other mandates isolating soft gluons from the perturbative 
kernel. At the same time, introducing the cutoff at any reasonable factorization scale 
eliminates logs as a large parameter in actual B decays: with maximal E 1 = rrib/2 
below 2.5 GeV the parameter ln-^ 1 can barely reach even unity. Therefore, we should 
not expect significant uncertainties from possible higher-order corrections. 

The basic features of the perturbatively generated distribution become manifest 
already in the model for the leading-log spectrum improved to include running of 
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a s , analyzed in Ref. 0. One has[] 

S ' A ' S(E y ) =e~ w ^ E -< ) (5) 



T d£ 7 d£ 7 

where 5 is the square of the Sudakov formfactor, and w{E 1 ) is the perturbative 
probability for a sufficiently hard gluon to be emitted so that the photon is left with 
energy below E 1 . Using the fact that soft gluons couple with a s {k±) and cutting 
the integration at k±<fi where /i represents the separation scale, the integrals are 
easily calculated 0: 

/r, x 8 /2tt 

■_in *ifa _i i„ «.w £ >„ 



1 J n My^fe) 1_ j n _ _9_ Y n a 

a 3 (m b ) a s (m b ) a s (» a s (y/em b ) 2w £ 



1-ln " 3 ^L, 



m b — 



a s (mi,J a s ( m W 27r — mi, 

(6) 

£ = m h — 2E^ . 

The representation of Eqs. ©-© has the added advantage that the integrated 
fraction $ 7 (-E) of the decay events with E 1 < E is directly given by 1 — S(E) = 
l — e~ w ( E \ The effective 'soft' coupling in Eq. (|) is known to two loops |§. 

The perturbative photon spectrum following from Eqs. was discussed in 

Ref. ||. The relevant point is that S(E) does not depend on /i as long as m^—IE > fi. 
This means that this part of the spectrum is shaped by gluons harder than /z; in 
practice this scale is around 0.9 GeV for E = 1.9 GeV. Therefore, the bulk of the 
perturbative spectrum below 2 GeV is indeed generated by sufficiently hard gluons, 
and can be trusted. Since the fraction $ 7 (-E) is still small, exponentiation of the soft 
emissions does not produce a radical numerical change. It is important, however, 
that this applies only as long as the cutoff in soft emissions is implemented. 

The upper part of the spectrum is mostly determined by the primordial distri- 
bution function F(x) with perturbative radiation affecting mainly the height of the 
spectrum and somewhat widening its shape. To a crude approximation, the central 
part of F(x) can be directly taken from the observed distribution at _E 7 >2.2GeV. 

The lower part of the spectrum is to some extent affected by the primordial 
Fermi motion as well, since <&v eTt (E) gets convoluted with F(x): 

$ 7 (£) = J dk + ^ cn (E-k + )F{k + ) . (7) 



lr This ansatz in fact produces even a somewhat harder spectrum than the one from the exact 
DL result based on soft factorization, which is evident when double emission is considered. If the 
soft gluons emitted individually lower E-y by amounts Si and 62 with probabilities dwi and dw2 
respectively, double emission would yield E 1 — ^—81—82 with the probability oc Aw\ &W2, while the 
ansatz places this contribution at E 1 equal to the lower of the two energies, E 1 = max{5i, #2}- 
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First, a long tail of F(x) to negative x can populate the domain of lower E 1 even 
in the absence of hard bremsstrahlung at all. This clearly would only enhance the 
estimate of $ 7 (2 GeV). Secondly, there is an effect of smearing due to typical k + ~ 
Aqcd- Since the perturbative spectrum decreases sharply with increase of ^ — E^, 
this normally enhances $<y(E) as well, at least with E in the deep perturbative 
domain - the perturbative distribution is much wider than the intrinsic primordial 
smearing in the 1 / uiq expansion. However, in practice typical momenta in F could 
be of similar scale as ^—E^; when k + kicks up the effective perturbative energy E 7 + 
k + /2 to the domain where $P ert is nearly saturated, the above approximation does 
not work anymore. Therefore, the primordial smearing can, in principle, decrease the 
purely perturbative estimate - in particular, when <&v eTt (E) is already significant and 
F(x) is a broad function. In practice, the nonperturbative smearing only enhances 
the tail of the spectrum in the relevant domain. 

The double log approximation turns out to be far too crude for B decays. This 
can be assessed by comparing the exact 0(a s ) spectrum [K| with the one obtained 
expanding the double log expression in a s : the resulting spectrum exceeds the com- 
plete expression by more than a factor of 3. Moreover, keeping double and single 
logs alone would yield an enhancement - rather than suppression - factor in the 
Sudakov exponent for E^ as large as ^y-— 75MeV. In this situation one better relies 
on the exact one loop spectrum: 

1 drP ert _ a s (2x 2 -3a;-6)a; + 2(x 2 -3)ln(l-s) _ 2£ 7 

r d£ 7 ~~ 3(i-z) ' x ~~m' ^' 

It is clear, however, that using a fixed coupling a s {m b ) would significantly under- 
estimate the actual spectrum. The BLM part of the second-order corrections has 



been computed in Ref. jTT| and was indeed found to increase the spectrum by about 
60% for low enough E 1 . 

There are subtle physical reasons for this BLM approximation to be of limited 
validity in the domain of relevant E T The BLM corrections for a given photon 
energy are obtained integrating over all gluon configurations consistent with the 
invariant mass of the strange quark-gluon system. In particular, the gluon can be 
energetic and collinear, or relatively soft yet emitted at large angles: 

M 2 X = 2m b (m b -E 1 ) « (9) 

The scale of the strong coupling, however is driven mainly by k±. Hence, the ef- 
fect of soft gluons is underestimated when averaged with those from short-scale 
configurations. 

This becomes important when double gluon emissions are considered. Their 
probability is proportional to a 2 s) and they are missed by the BLM terms. In contrast 
to ordinary corrections, this effect has an extra enhancement factor of ln(l — x)/x 
for each emission. Thus the common wisdom about the BLM dominance is not 
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expected to apply here. And there is an additional numerical enhancement: dou- 
ble (or any multiple) emission allows lowering E 1 by the same amount with softer 
gluons (e.g. lower k± for same uj); since a s grows fast, this would yield a significant 
increase of the radiative tail overlooked by the naive BLM scale fixing. This can 
also be understood in a complimentary way applying the OPE to soft perturbative 
bremsstrahlung. 

This problem obviously does not arise if the Wilsonian prescription for separating 
long and short distance dynamics is implemented when computing the perturbative 
spectrum to avoid double counting of soft gluon modes. Then the infrared domain is 
excluded and the coupling never grows too large. Moreover, soft gluons are simply 
excluded from computations regardless of the strength of the coupling. This elimi- 
nates both the possible enhancement in the soft emission amplitudes and strongly 
suppresses multiple emissions altogether; for the latter would now lower E 1 by too 
large an amount. This approach allows using a single-gluon-generated spectrum 
and makes the BLM-type accounting for the running of a s trustworthy. 

The estimated fraction $>v eTt (E) depends mainly on the exact value of the b quark 
mass. Actually $P ert is primarily a function of ^ — E rather than of E itself. It is 
evident that one should use the perturbative running low-scale mass m fe (/z) here, not 
the pole mass: excluding emissions of soft gluons requires - due to gauge invariance 
- discarding analogous virtual corrections, including the soft corrections in the b 
quark self-energy. Then the pole mass never appears in the problem, and the heavy 
quark is seen in perturbation theory as a particle with the mass near m^/i). 

In principle, the Wilsonian prescription can be introduced in different ways. A 
consistent scheme suitable also for genuinely Minkowskian processes and beyond 
one loop has been elaborated |TB[ . In essentially one-loop corrections (including 
arbitrary order in the BLM resummation) it reduces to a cutoff in the gluon energy 
uj > fi (in the perturbative computations) in the heavy quark restframe. The value 
of the "kinetic" running quark mass and its running is well known (for a review see 
Ref. 1): 

m fe (l GeV) = (4.57 ±0.05) GeV (10) 

(a somewhat larger value - 4.63 GeV - however with twice the uncertainty was 
reported by CLEO 0). A subtlety should be kept in mind though: with such a 
cutoff scheme there are different short-distance masses perturbatively related to each 
other. The one entering here is m (/i), the mass defining the rest energy, rather than 
the kinetic mass in Eq. ([T0|). The former slightly exceeds m;,(/z): 



kin, - , 4a s / , , (3 Q a s 



m (/i) = m kin (/i) + -— fi 1 + 



9 7T V 2vr 



m b 13 

m 

2/x 6 



+ -) +o M?)- < u » 



The difference constitutes about 60 MeV at /i = 1 GeV. (The straightforward all- 
order BLM summation yields a larger shift, however for consistency we retain only 



terms through (3qoi 2 s .) With such the ratio ^ 1 is about 0.86 at E~ = 2 GeV. 



nib 



The perturbative expressions existing in the literature [0, 1 1] for the one-loop 
spectrum do not include the effect of an infrared cut, even though this computation 
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is not difficult. For low enough E 1 this would not affect the single-gluon spectrum 
as have been already noted for the double log corrections. Strictly speaking, a gluon 
with u<n can bring as low as //, so the cut affects the spectrum in a wider 
domain. This is possible, however, only when the gluon and the s quark fly back-to- 
back, a rather rare and - being governed by a short- distance coupling - suppressed 
configuration. The enhanced kinematic sensitivity is absent as soon as the gluon 
is emitted into the forward hemisphere. Hence, to a good approximation the effect 
of the cutoff can be neglected at E 1 below mb ~ M as it happened in the double log 
ansatz flU). 

Above this borderline the perturbative spectrum is modified more significantly, 

2 

its growth slowed down with the spectrum vanishing around ^ — It also has 

the two-body 5-function at E 1 = m ^ suppressed by a finite amount by virtual 
corrections. 




E, GeV 



1.5 1.6 1.7 1.8 1.9 2.1 

Figure 2: Perturbative fraction of decay events $P eTt (E) having E 1 <E. 

The estimated fraction of the decay events generated solely by perturbative glu- 
ons is shown in Fig. 2. We see that one expects $P ert (2 GeV) to be about 6%. 

2.2 Combined spectrum 

The overall spectrum is obtained by convoluting the perturbative contribution with 
the light-cone distribution function E(x). Hence it requires the perturbative part in 
the whole domain including the interval irib — 2/x < 2E 1 < rrib where it is modified by 
the gluon cutoff fi. However, the exact details in this relatively narrow domain are 
not too significant due to primordial smearing, as long as the overall normalization, 
$F> crt (-E) = i a t E > ^ is maintained. Using the double log computation as a guide, 
we extrapolate <&v eTt (E) linearly above E = mb ~^ ■ the end-point 5-function in the 
spectrum is fixed by the normalization constraint. This approximation is expected to 
slightly underestimate $P ert in the lower end of the window ^ — fi, and overestimate 
it near the maximal E 1 except the endpoint itself. The bias in the resulting total 
fraction should then be insignificant and presumably only underestimates $ 7 (£') at 
E near 2 GeV. 
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As illustrated in Ref. 0, while the purely perturbative spectrum and $P ert (£' 7 ; /x) 
depend strongly on the choice of /i, the overall $ 7 (-E) comes out fairly insensitive to 
it. For the ji dependence of $P ert is strongly offset by the renormalization of mj(/i), 
i.e. by the shift in the starting point of the parton-level spectrum. Therefore we 
can infer numerical estimates using a particular choice of the cutoff scale, namely 
/i = 1 GeV for which the nonperturbative parameters are routinely fixed. 
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Figure 3: Fraction of decay decay events <J? 7 (2£) including both perturbative 
bremsstrahlung and primordial Fermi motion. Solid and dashed lines correspond to 
ansatze Eqs. (|T^ ) and ([i~4| ) for the light-cone distribution function, respectively. 



The expected fraction $ 7 (i£) is shown in Fig. 3 and the spectrum itself in Fig. 4; 
we superimposed the experimental points onto our predictions in the latter. We 
took the ansatz of Ref. [H for the distribution function F(x), but simplified it to an 
essentially a single-parameter one: 

F(k+) = N (k + -K) a e~ ck+ 9{k + -K) . (12) 

c gauges the scale of the light-cone momentum, and a determines the dimensionless 
ratio ^7 = Numerically a = 2 yields a value very close to the experimental 

ones with /4(lGeV) = 0.43 GeV 2 and A(l GeV) = 0.65 GeV. It is curious to note 



Pi 



for arbitrary a. Such a value is favored by theory and would hold exactly if the SV 
sum rules were saturated by the P-wave states at a single mass. 

To check the sensitivity to the shape of F we have also employed the alternative 
ansatz 

F(k + ) = N(k + -Afe- {d{k +- R))2 0(Jfe+-A) (14) 

with the tail fading out much faster, and (3 adjusted as to reproduce the same kinetic 
expectation value for a given m^. (At a=2 this corresponds to (3 ~ 0.658 and yields a 



that ansatz ([12] ) yields for the Darwin expectation value 

(13) 
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Figure 4: Expected photon spectrum for the two ansatze for F(k+). The relative nor- 
malization of experimental points and theoretical predictions is not fixed. 



factor of 2/3 smaller Darwin expectation value). As expected, is practically 

insensitive to the concrete ansatz; even the difference in the spectrum itself lies 
below the effects of higher-order nonperturbative corrections. It should be noted 
that the predictions shown in Fig. 4 do not include Doppler smearing due to the 
initial velocity of B mesons produced at T(4S f ). Assuming absence of a correlation 
of experimental measurements with the direction of the decay relative to the B 
momentum, the additional broadening is easily incorporated. The corresponding 
plots are given in Fig. 5. 

Therefore, we arrive at the conclusion that about 12% of decay events should 
have E^<2 GeV. Half of this probability is generated perturbatively, and it already 
exceeds the 5% of CLEO's fit. The fact that CLEO's fit represents too narrow a spec- 
trum can be inferred also from the size of the second moment or dispersion through 
an alternative perspective: the value quoted by CLEO is essentially oversaturated 
by perturbative contributions alone. 

It should be noted that although the fit was used in the analysis of Ref. [Q, 
the total decay probability was obtained there adopting $ 7 (2 GeV) = 8.5t|f% from 



Ref. []I4"| , even though the fit itself yielded about 5%. 

At present the part of the spectrum below 2 GeV does not seem very certain. 
If, however, future measurements establish that significantly less than 10% of the 
spectrum resides there, this would send the alarming message that the standard 
treatment of strong interaction effects - in particular the gluon corrections as they 
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1.5 2.5 3.5 4.5 

E r (GeV) 

Figure 5: Photon spectrum in the lab frame where T(4S I ) are at rest, including Doppler 
smearing with /3b = 0.067. The relative normalization is arbitrary. 

are routinely applied to decay distributions for beauty hadrons - is unreliable. This 
in turn would cast doubts on the numerical accuracy of extracting m&, etc. from 
the photon spectrum and even more on relating the lepton endpoint spectrum in 
b — > utv to the b — > s+7 spectrum. 

It can well be that the fully integrated spectrum is still robust against such 
effects; however the higher moments and the shape in general would be affected. If 
CLEO's fit were to be vindicated by future data, one could not extract the quantity 
A (let alone the kinetic expectation value) using the quoted formulae; for those 
expressions were based on the assumption that gluon radiation can be described 
perturbatively over the whole range up to _E 7 = 2 GeV. 

The simplest resolution of a possible discrepancy between the CLEO spectrum 
and the theoretical expectation is that in particular the data point at the 1.9 GeV 
bin represents a low fluctuation. Future data will settle this issue. 

Of course, the first and, in particular, second moment of the photon energy 
decrease if a spectrum closer to QCD expectations is used. Unfortunately, one 
cannot re-evaluate them based only on the data points given in Ref . || . Extracting 
the value of the b quark mass it relied only on the moments evaluated over the range 
E 1 >2 GeV, hence there was a limited sensitivity of the purely experimental input 
to the lower part of the spectrum. However, this simply relegates the sensitivity 
to actual physics of gluon bremsstrahlung to theoretical evaluation of such defined 
moments. 
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3 Skeptical remarks 



One of the practical applications of measuring the prompt photon spectrum is an 
attempt to estimate more reliably the fraction of semileptonic b — > u decays expected 
theoretically with Ei in the end-point region beyond the kinematic range for b — > c. 
While appealing as a nontrivial physical statement, in practice the relation between 
the integrated semileptonic end-point spectrum 



T s \(b — > u) Je dE, 
and the corresponding weighted B — > X s +7 spectrum 



1 ,„ (M B \ dr 6s7 



v b ,.,J E AE A~f- E ^ (16) 

is rather fragile - a fact emphasized already in the first dedicated paper |T3: the 



actual value of m;, turns out to be insufficiently large to make higher order cor- 
rections irrelevant, and the allowed interval of lepton energies Ei > | M D is too 
narrow to suppress higher-order corrections violating universality of the Fermi mo- 
tion. There are two basic nonperturbative effects which are potentially dangerous 
here; we illustrate them in more detail than it was done in Ref. [|T^] . 

The first one is related to a particular class of l/m| corrections associated 
with four-fermion operators responsible for generalized Weak Annihilation (WA) 
effects. While typically at a few percent level in fully integrated B widths, they 
are concentrated in the end-point domain in semileptonic decays and thus are en- 



hanced by an order of magnitude for the narrow slice in Ei fL6fl . Their size is 
controlled by the poorly known nonfactorizable terms in the expectation values 
{B\fry^{\ — 7 5 )mr)v(i — 7 5 )6|£>); they are different in charged and neutral B and in 
general do not vanish even in B°: 

XV (h \~v*(h ^ 327^ 2 {B^l-^uu^jl-^B) 



The nonfactorizable contributions were addressed in Ref. [17]; the non- valence oper- 



ators, however, are rather uncertain. The reasonable estimate of Ref. |T8[ suggests 



that this effect constitutes 10 4- 30% of the semileptonic width with Ei > 2.3 GeV. 



A more detailed analysis can be found in the original papers |16|, |TI 

No such effect arises for the B — > X s + 7 width at tree level. It is generated to 
order a s ]15| and is given by a different four-fermion operator. The correction to 
T(B — ► X + 7) is expected smaller here; a more careful consideration reveals that 
the end-point fraction can be affected, yet less significantly than in B — -> X u tv. 

The second problem is of a general origin and related to the fact that the interval 
in the electron (or photon) energy Ei > 2.3 GeV is too narrow in practice |15||. To 
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quantify this observation, let us recall the I/tuq corrections in the b 
spectrum |19|, [3J : 



u lepton 



J_dT k = 2 /2£,f 
T dE e \m b J 



4E, 



m b 



10E P 



mr 



3m, 



Ej 



3m\ 



(m b -E e )5'(^-E e 



a4+ 



+ 



10E e 
3mf 
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12m b 



and examine the effect of the chromomagnetic operator. It has a negative 5-function 
contribution in the end point: 



i dr sl 
r dE t 



2 1 




+ term oc fi w 



(19) 



The spectrum can be neither singular nor negative, of course. What this means 
is that one has to average the lepton energy at the very least over an interval 
Amin = j2 ^ — 70 MeV; a similar shift and a need for smearing applies a priori to 
b — > s + 7 events as well. On the other hand it is quite evident that this effect is 
completely missing from the Fermi motion and from the distribution function F(x) 
- in contrast to the latter it is a non-universal pro cess- dependent nonperturbative 
correction and would even have the opposite sign in decays of B*^\ Thus it is absent 
from the naive leading in l/m b relations between the observables in Eqs. ( |T5"1) and 



(0) 



The shift by 70 MeV does not affect much the total decay rate - yet it is very 
significant for the narrow end-point slice: for the parton spectrum it would decrease 

This is more than a 50 percent reduction for 

E- 



the rate by a factor 1 — if — , Mg ocn 

: 2.3GeV even with conservative assumptions about m b . 



We can suggest a simplified way to improve the end-point relation between the 
photon and lepton spectrum close in spirit to the proposal to be made below in 
Sect. 5.2 for the invariant mass M\ distribution. To account for the chromomagnetic 
interaction effect in Eq. (plj) one can make a shift in the lepton vs. photon energy 

by a formally 0(l/m b ) amount 8Ei 



2Mg. 
3 mi, ' 



r s i(6 



u 



dE, 



u 



dE, 



oc 



1 , M B 

1 12 



dE^ 



M 



B 



E, 



dT 



bs'y 



dE^ 



(20) 



Physically this is motivated by the different interaction of fast light-like light quarks 
with the chromomagnetic field in the two processes leading to a relative energy shift. 
However this is not a rigorous prescription and its accuracy is difficult to quantify. 



2 Chromomagnetic field B(0) emerges as a commutator of different spacelike momentum op- 
erators of b quark, while Fermi motion depends only on a single space momentum and energy 
operator. The term oc /i 2 is the first term in the effect of the primordial motion. 
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Since the effect is significant, it may introduce an essential uncertainty even with 
this improvement. 

It is important to keep in mind that the above constraint on the minimal resolu- 
tion is of a rather fundamental nature which goes beyond a conceivable improvement 
of the Fermi motion description M . Knowledge of the subleading distribution func- 
tions would not enable one to shrink this interval. Of course, the literal expression 
for higher-order terms like Eqs. (|T7|), ( |iT?D or (^) going beyond the leading- twist 
expansion in this end-point domain cannot be used literally to improve the naive 
leading-order expressions, since they are superimposed onto the nonperturbative 
leading-twist Fermi motion. In principle, it cannot be excluded that the latter 
somewhat softens the numerical instability of the end-point expansion - yet one 
cannot count on a radical improvement of the accuracy: the fact that allowing only 
a 70 MeV interval of the end point energy yields a 100% error in the spectrum in 
the best scenario, is a benchmark which is difficult to get around. 

From the general perspective the problem is related to presence of another pa- 

2 

rameter in the problem - the ratio of the scale — to the typical hadronic mass. 

„2 

With this scale amounting to only a few times ^ the relation for the end-point 
spectrum is violated by a fixed amount even for asymptotic values of m&. The ac- 
curacy of the relation is not governed solely by l/m b as is usually stated, and the 
actual numbers favor a possibility for significant corrections. This point is often 
missed in the literature. 



4 Improving the M\ distribution in B — ► X u tv 

The integrated width T(B — > X u iv) is under good theoretical control |20|1 . Yet to 



measure it one has to disentangle it from the dominant B — > X c lis width. Vertex 
detection is not so highly efficient that one could achieve this goal solely or even 
mainly by rejecting the b — > c background through finding the secondary charm 
decays. One has to apply cuts of a mainly kinematic nature. The most direct 
way is to consider the hadronic recoil mass spectrum and impose a cut in the 
invariant hadronic mass Mx <Md, which has the advantage of retaining almost all 
b^u events. For experimental reasons one will have to lower the cut to confidently 
reject the b— >c background. The situation is quite favorable since even a cut as low 
as 1.5 GeV leaves more than half of the signal events. Yet one still wants to place 
it as high as possible, preferably above 1.65 or even 1.7 GeV. Considering the huge 
statistics that can be accumulated at modern B factories the central issue becomes 
how well one can theoretically evaluate the fraction of the b —>■ u width falling into 
a restricted domain of phase space. Raising the cut reduces the rejected fraction 
and thus suppresses the theoretical uncertainty in the measured rate even without 
precise control over the fraction. 

While the M^-spectrum is sensitive to the values of the nonperturbative pa- 
rameters as well as the heavy quark distribution function F(x) describing Fermi 
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motion, the dependence on the latter is moderate if the current knowledge and the 
constraints on F(x) are incorporated in full. This particularly applies to the par- 
tially integrated spectrum. Following Refs. [2T| we consider the fraction of events 
with the hadronic recoil mass below M max : 

M M) = -^I\m x «£^ (2!) 
r s i(o — > u) Jo dM x 

Even without a dedicated analysis of all the uncertainties it is evident that $ s i(1.7 GeV) 
must exceed 0.55 and be below 0.9. This observable is thus known with at least ±24% 
accuracy translating into a ±12% error bar in \V u i\. Even allowing for additional 
stretching of uncertainties this would represent at the moment a good benchmark 
for quantities relevant for \V u b\. We found no substance to claims of unrecoverably 
large uncertainties in this observable. 

Ultimately we would like to be able to decrease the theoretical uncertainties 
in $ s i(M) further still - and achieve it without having to push the cutoff M too 
close to Md or degrading the confidence in the assessed error bars. One obvious 
possibility is to introduce an additional cut rejecting small invariant mass q 2 (q 
denotes the momentum of the lepton pair) where the strong interaction effects on 
M\ are maximal due to the significant recoil.^ Choosing, for instance q 2 > 0.35m 2 . 
practically eliminates the impact of the primordial Fermi motion on Mx < 1.7 GeV 
events. However, imposing additional cuts has the serious drawback of introducing 
a significant dependence on higher-order effects which are difficult to control and 
which are all too often missed in estimates. The overall energy scale governing 
the intrinsic 'hardness' of the decay rate deteriorates being driven at large q 2 by 
m b — Vq 1 instead of m& itself [13| . Fortunately, there is an alternative way not 
employing additional cuts. 

The idea is based on the realization that the dangerous small-g 2 domain in B — ► 
X u lv decays is nearly identical to the case of B^X S +^ in respect to the hadronic 
mass distribution, at least in the limit of large ray,. Moreover, its accuracy is governed 
by a genuinely l/my, expansion. Therefore, we can directly utilize experimental 
information about the B— >X s + / y spectrum to control this kinematics, in a far more 
reliable way than the charged lepton end-point spectrum. 



4.1 Scaling behavior and universality relations 

To illustrate the idea, let us consider first B^X u £u at q 2 = 0. The concrete Lorentz 
structure of the decay vertex {V — A in b-^utv and tensor in b— > s + j) does not 
matter in the heavy quark limit; the effects of strong interactions rather depend on 
kinematics which are equal for vanishing invariant mass of the particles recoiling 
against the final state hadronic system. The decay distributions are governed by 

3 A cut, say on the lepton pair space momentum \q\ achieves a similar purpose since it is 
kinematically related to M\ and q 2 . 
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the light quark propagator emitted with light-like momentum \k u \ ~ ^ f 1 — ^ 
Then one has 

1 dr sl (g 2 = 0) 2 _L_ ^dTV / Mf-ivr| \ , , 

r sl (g 2 = 0) dMi {Mx> 2M B T bs , dE, { ™* J ' ^ ] 

since kinematically E 7 — M - B 2J ^ adr . Here and in what follows we assume the b^u 
transitions when refer to semileptonic decays. Thus, at least for q 2 = we can 
directly obtain 1 — $ s i(M) by integrating the photon spectrum in B — > X s +7 up to 

L- = M|-M2 

2 



Relation does not hold when g in the two processes are different, namely 
away from the point q 2 = in the semileptonic decays. However, there is a scaling 

over the variable M^ adr / ^1 — reflecting the universality of the Fermi motion 

15fl . It simply states that the hadronic mass squared is proportional to the decay 
light quark momentum \k u \: 

_L dTsi^g 2 ) , 2 , _ 1 m 2 , 1 djj^ / Ms M | 



T sl (q 2 ) dM 2 1 x> 2M B m 2 -q 2 T bs ^ dE^ \ 2 2Mb J ' 1 ' 

This relation deteriorates for larger q 2 and breaks down when m 2 — q 2 becomes too 
low. However, at q 2 constituting a significant fraction of M B the effect of Fermi 
motion on the integrated fraction $(M) disappears which simply becomes unity in 
this approximation. 

Having this in mind, we can implement the scaling relation in practice to evaluate 
$ s i(M). For the required rate with Mx exceeding M we have 



/ - > 1 dr si f llf2 1 d 2 r sl (M 2 ;g 2 ) 

with the last differential width in the scaling approximation given by Eq. 
^T s i(g 2 ) and ^t 1 denote the same total fixed-g 2 width.) Then we have 



M 



B ml-q 2 

M l- M \ , , ,, i dr., t 2M.' i dr 



i-*, l( M; X )=/ - deW)^-/ 2M " ds,=i-^, (25) 

r s i dg^ Jo r fcS7 d^; 7 



where x(<? 2 ) is a weight function one may (but does not have to) introduce to suppress 
contributions from small g 2 kinematics. Straightforward arithmetics yield 

m b m 2 

1 ~ $ si(M; X ) = I 2 2Ms dE, 0(£ 7 , M; X ) J- j (26) 



with 



„ M 2 mJ 
2 0. 



0(£ 7 ,M; X )= / 6 I7 ^ i ^dg 2 X (^ " 
Jo 



/o * ^ ' T s i dg 2 
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For not too large q 2 relevant in Eq. ( p5| ) the parton approximation 



r m ^ = 6 fl-il) 2 -4(l-ilV (27) 
T sl dq 2 V m£/ V m% 



is quite accurate (see Ref. pjj). In what follows we assume xiQ 2 ) = 1 ( no additional 
discrimination at all); then we arrive at the following universality relation 



d> sl (M) = / 2 2Mb dE y 0(£ 7 , M)— — _£i£2 ( (28) 

JO 1 b S7 



d£ 7 



2r d r 4 2.EL M 2 

«^ = 1 -(i=^ + (i=^- ^ = Mb ' r = ss- 

It allows one - as promised - to express <3> s i(M) directly via the photon spectrum. 
The weight functions 4>{E^) depend on the chosen M^-cutoff M and vary from 

1 at <C Mb / 2 down to zero at the maximal E 1 = ^ M where the single point 
q 2 = contributes. They are shown in Fig. 6 for M = 1.6 GeV, M = 1.7 GeV and 
M= 1.8 GeV. As anticipated, at E 1 < 2 GeV the weight <p(E^) is practically unity for 
all reasonable cut masses M. This means that the low part of the photon spectrum 
contributes only as a total integral, in other words, enters as the overall fraction of 
events discussed in the previous sections. The spectrum at E^ > 2 GeV enters with 
decreasing weight: 

M B M 2 

1 - $ sl (M) ~ $ 7 (2 GeV) + / 2 2Mb dE y <j>(E,; M)— — . (29) 

J2GeV 1 bs-y dE^ 

The practical importance of clarifying the actual value of $ 7 (2GeV) becomes then 
manifest . 

The consequences of the Fermi motion description of the large-recoil light-cone 
kinematics similar to explicit relations (E3) have been discussed in the literature more 



than once (see, e.g. the recent paper p2|). Taking a different perspective, though, 



it is usually claimed that the outcome for $ s i suffers from significant uncertainties. 
We do not share this pessimistic viewpoint for actual B decays. For example, the 
obvious constraint that $ s i(M) is exactly unity regardless of dynamics at sufficiently 
high M approaching Mb is usually neglected. As illustrated earlier in this section, 
in reality we need to evaluate only the ah initio suppressed fraction 1 — $ s i(M), and 
therefore even its crude estimate does not introduce a significant uncertainty in \ V u b\. 

Universality relations and the scaling relations in general are clearly not exact. In 
the given form they are robust against perturbative effects; in particular the Sudakov 
double log corrections would not violate them. Yet they hold only to leading order 
in l/m&. For example, terms of order A were neglected in M\ compared to the 
leading ones ~ m^A. Since 1 — $ s i(M) itself is a small fraction and thus introduces 
only a correction in extracting \V u b\, its approximate knowledge already suffices for 
our purpose. Nevertheless further improvement can be achieved by evaluating l/m& 
corrections to the universality relations Eqs. (E 
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1.6 1.8 2 2.2 2.4 

Figure 6: Weight functions (/)(Ej) for three values of the cut on hadronic invariant mass, 
M = 1.6GeV (dashed), M = 1.7GeV (solid) and M=1.8GeV (dot-dashed). 

4.2 Leading l/m^ corrections 

Accounting for the A terms referred to above is quite simple, and in general they 
can easily be incorporated in both semileptonic and radiative decays. However, this 
would not give the full answer. The universality relations make use of two basic 
facts, namely the process-independence of the leading-twist Fermi motion in the 
light-cone kinematics and the scaling of M\ (or a certain function of it) with 1— 

While it is easy to account for the l/m b corrections to the former, l/m b scaling 
violation is governed by independent distribution functions. 

There are two circumstances which help to overcome these obstacles. In practice, 
only a limited range of q 2 yields the unwanted background due to Fermi motion; 
therefore the scaling violation in passing to q 2 > is not dramatic and can be 
dealt with. Besides, we do not aim at detailed distributions, but rather need an 
integrated probability. In addition, in the small-g 2 domain where the Fermi motion 
is potentially most important, its universality is ensured by the full b quark mass 
rather than a lower momentum scale. Below a way is suggested to implement these 
ideas in practice. 

Instead of assuming the leading twist scaling in M\j{\ — for decay distribu- 

b 

tions, we can employ its l/m b improved ansatz: 

_i dT sl (fi 2 ;q 2 ) 1 m 2 - q 2 1 &T bsi / Mb ^ m 'j A ( . 

T sl (q 2 ) d/i 2 " 2M B m 2 T bsj d£ 7 I 2 2M * ^F?> ) ^> 

with fj, 2 given by 

" 2 =^+^) + K*) Mi+c (*)3 (3i » 

where A ~ Aq CD , B ~ and C ~ 0(1) are effective parameters accounting 
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for l/m b corrections. In other words, we assume fi 2 instead of M\ itself to scale 

2 

proportional to 1 — fy. 

b 

In the simplest parton model one has, as an example, 

a = r2 > B = jr 1 -r^> c = 77^-2-> x = -n*> ( 32 ) 

M B 1 — x (1— x) 2 M B 

which would incorporate kinematic effects. This is analogous to deep inelastic 
lepton-nucleon scattering, where the introduction of Nachtmann variable leads to an 
extension of the scaling regime. Simultaneously, this can also account for corrections 
to the heavy quark symmetry for different weak decay currents (V — A in b — > uiv 
vs. tensor in b — > s+7). We do not rely on a model here, but rather can fix the 
parameters theoretically minimizing the strong interaction corrections computed at 
given q 2 in the standard l/m& expansion JlJj, |T|. This determines the ansatz pa- 



rameters in a model- independent way. After that, the straightforward integration 
would yield the l/m& and perturbatively improved version of the universality 
relation (p8|) . 

A clarifying note is in order here. The standard l/mb expansion yields non- 
perturbative effects only as an expansion around the free quark kinematics, e.g. 
M\ = (Mb — mb) 2 + (Mb — m^rribil — for b^>utv consisting of the (^-function 
and its higher derivatives. Their direct integration in the observables of interest 
makes little sense in the present context (likewise, it cannot be used to quantify the 
accuracy of various constrained inclusive probabilities in the case of limited energy 
release, as can be sometimes seen in the literature). Instead, here we can compute 
at given q 2 the hadronic moments ((Af|), (M^), • • • ) in terms of heavy quark pa- 
rameters uif,, a4> ^g; Pdj ■ ■ ■ directly and in the ansatz (j30|)-(|3ll), and require them 
to match. In this way 1 / corrections are properly accounted for. 



5 Conclusions and outlook 

By the turn of the third Millennium \V u b\ had been determined through a measure- 
ment of the inclusive rate T(B — > X u lv). Most recently - after this paper had 
largely been completed - a new result from CLEO has appeared with an ex- 
traction of \V u b\ from the end point lepton spectrum. Yet the theoretical reliability 
level of these methods varies considerably. The inclusive rate T(B — >X u iv) suffers 
from little model dependence. The main problem arises in distinguishing it against 
the dominant b^c background and how to model the latter 

In our considered judgement measuring the hadronic recoil mass spectrum in B — » 
X u lv will provide one with the most reliable extraction. Placing a high premium 
on model-independence we think that a measurement of the lepton energy endpoint 
spectrum per se cannot provide a competitive extraction of \V u b\ as long as a cut 
around Ei ~ M ^ Af M ° remains, whether or not the B — > X s + r y spectrum is known. 
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Nevertheless studying the end point lepton spectrum is of high theoretical in- 
terest. The fact that the analysis of Ref. [0] yielded a value for \V u b\ very close to 
other determinations places bounds on the expectation values of four-quark opera- 
tors. Naive estimates result in a typical limit around 



(B\H(l- l5 )u Ull (l- l5 )b\B) 



2M 



B 



< 0.012 GeV 3 



(33) 



(an average over charged and neutral B is implied) which lies in the range of previous 
estimates |1| [L7|, |18| . A more careful analysis is required, though, to obtain 
reliable numbers. In particular, an improvement in treating the end-point relation 
between the semileptonic and radiative spectra may turn the bound into a nontrivial 
measurement of these expectation values. The interpretation would be much better 
supported if measurements were performed separately for charged and neutral B 
mesons. 

High statistics studies of the hadronic recoil mass Mx in B — > X tv provide 
good means of obtaining an accurate value for \V u b\. The photon spectrum in B —> 
Xs+j decay gives a more or less direct measurement of the heavy quark distribution 
function, far beyond indirect constraints on its gross features inferred from estimates 
of a few moments available a few years ago. The universality relations allow to 
convert the spectrum into quantitative estimates of the required fraction $ s i of 
B — > X u tv decays with high M x rejected to suppress the b — ■> c background. For 
example, using the simplest version - Eqs. (|28|), (|29|) - we obtain $ sl (1.7 GeV) ~ 0.65 
assuming $ 7 (2 GeV) = 12%. We have argued that already the present uncertainty in 
this b— >u fraction is not significant, and suggested the way to improve it for further 
studies by incorporating l/m& effects. 

There are statements in the literature that extracting \V u b\ 
bution is seriously plagued by uncertainties in the heavy quark distribution function 
F(x), while decay rates with a cut q 2 > (M B — M D ) 2 are accurately calculable and 
stable against higher order effects. Such arguments miss, however, some basics of 
the heavy quark expansion and in particular how the Fermi motion itself emerges 
in the OPE. The claim that some fraction of T s \(b — > u) constrained by a cut on 
q 2 does not depend on F(x) is actually justified only in a negative way: this width 
is governed by a whole series of unknown expectation values of higher-order heavy 
quark operators which cannot be resummed or otherwise related to F(x). The lat- 
ter is known reasonably well since it is strongly constrained by sum rules and can 
directly be bounded by data. The corrections to widths with a cut in q 2 are obscure 



from the M\ distri- 



in this respect; there is actually indirect evidence that they are quite significant |L8 
A similar reservation was expressed recently in Ref. |26[ . 



To give a brief illustration of missing pieces, let us consider a constrained fraction 
of the B— >X s + / y events 



l-$ 7 (£) 



1u b i-i-, 

~ , „ dr 



dE. 



bs'y 



(34) 
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Ignoring the perturbative corrections, the usual l/m& expansion always yields the 
spectrum in the form of expanding around the free-quark kinematics: 



1 dTbs-f 



aS^-^) + b5'{E 1 -?f) + C ,T(£ 7 -^) + ... , (35) 



where a, b, ... are given by the expectation values of local 6-quark operators over 
B. Naively computing 1 — $ 7 (i?) in this way would yield unity for any E> ^ - the 
result clearly unjustified until 'hardness' Mb — 2E is high enough in hadronic scale. 
Moreover, if $ s i(M) were computed in the same naive l/m& expansion as was used 
to quantify the corrections to the g 2 -constrained semileptonic widths, one would 
have $ s i(M) = 1 for any M ^ 1.8 GeV, with the decrease only due to perturbative 
effects. 



Similar reservations apply to the computation [IT] of 'incomplete' moments eval- 
uated over the restricted domain E 1 > 2 GeV which have been used by CLEO to 
extract m&. The offered expressions depend on the value of the lower cutoff (2 GeV) 
only in the perturbative piece. Switching off the perturbative corrections, there 
would be no dependence on the cutoff, which cannot be true in general and thus 
introduces an additional uncertainty. It is likewise obvious that the true dependence 
is more complicated and involves both perturbative and nonperturbative effects in 
a more intricate way than a simple sum. This complication would be reduced if the 
immediate part of the spectrum below 2 GeV can be incorporated into the analysis. 



In radiative B decays one indeed expects only a small fraction of the spectrum 
to reside below 2 GeV, yet not an insignificant amount: we anticipate about 12%. It 
contributes a barely significant amount to the overall BR(5— > X s + / y). Nevertheless 
its accurate measurement is important for various reasons. We expect - and hope - 
that new data sets with even higher statistics will reveal a low energy spectrum at 
the predicted level. Detailed study of this part of the spectrum are instrumental for a 
comprehensive program of exploring local quark-hadron duality; here it is probed in 
the somewhat special case of hard heavy-to-light flavor transitions which has its own 
advantages as well as complications. Should it be confirmed, we expect the (properly 
defined) heavy quark parameters A(l GeV), fi^ (1 GeV) - and possibly p 3 D , p\ s - to 
be measured with the comparable precision and to add confidence in our present 
knowledge of these parameters. Moreover, combining the whole set of emerging 
new data with tight theoretical constraints and precisely extracted mb(lGeV) will 
allow one to address in a quantitative fashion higher-order corrections - a possibility 
considered so far less than remote in conventional approaches. 

One of the important quantities in this quest is just the integrated low-i^ fraction 
of radiative decays $ 7 (2GeV) - even leaving aside the aspect that it is directly 
related to the determination of $ s i(M) at experimentally relevant cuts M« 1.7 GeV. 
Finding in the data a fraction of only, say, about 5% would tell us that some basic 
elements of the description had to be revisited. A possible interpretation would 
be that processes involving the emission (or excitation) of "real" gluonic degrees 
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of freedom start to exhibit perturbative local duality at higher scales than when 
only quark degrees of freedom are involved. This actually represents a conceivable 
scenario [J2U|J . The higher onset of duality here would mean that a perturbative 
treatment of gluon radiation is not reliable till one reaches a higher scale. From 
the specific case under study here with the separation scale around E 1 ~ 2 GeV one 
would be sensitive to the duality onset scale near 1.5 GeV; the lower part of the 
spectrum would allow to probe duality for even harder gluons. The total integrated 
decay rate still can be fairly insensitive to such effects. The photon decay spectrum 
from this perspective is a complementary tool to study details of local parton-hadron 
duality compared to hadronic mass distribution in semileptonic b — > c decays J20] . 
It is more ambiguous as a probe of lower scales due to more significant interference 
with the nonperturbative Fermi motion, yet may improve at lower E 7 if prompt 
photons and their production mechanism can be reliably identified. 

Although physics of gluon bremsstrahlung is quite different in semileptonic and 
radiative B decays, finding a significant violation of gluon-hadron duality in the 
latter at E 1 <^ 2 GeV would raise concerns in treating, say the hadronic moments 
in the former. This especially applies to higher moments in the environment of 
present CLEO measurements: the cut on the lepton energy Ei>1.5GeV effectively 
decreases the expansion mass parameter down from the usual energy release mj- m c . 
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